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In the 10 x 10 array on the cover, 22 of the cells 
have been selected by the following scheme. Two-digit 
random numbers were drawn (the numbers used for this 
example were taken from the table of random digits given 
on page 13 of PC33, two digits at a time, reading down 
the page). The first number was 57, which indicated the 
first cell of the cover pattern. Successive random 
2-digit numbers point to other cells, but a "hit" can 
only be scored in a cell adjacent to one previously hit. 
Thus, at the start, only cells 47, 56, 58, and 67 are 
eligible, and the second hit is scored when one of these 
(namely, 67) appears in the random stream (which was the 
14th random draw from the table being used). At that 
point, the number of eligible cells enlarges to include: 


His 5a 58, 66,68, and 77. 


This procedure continues until all the cells are 
hit; the score is then the number of random numbers used. 
For a 10 x 10 array, this is likely to take over 900 
numbers. To fill the 22 cells shown on the cover, it 
took 150 random selections of 2-digit numbers. 


For other size arrays, each game would take a 
different number of draws of the random numbers, and the 
result, for any given size of array, is a distribution. 
The accompanying table shows the results for 100 games 
played with arrays of size 2, 3, 4,...,20; the calculations [> 
were made by Associate Editor David Babcock. 
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The table shows that there is a wide range between 
the smallest number of draws needed to fill an array and 
the largest number. The average number of draws plots 
as a fairly smooth curve. When plotted on log/log 
paper, the points appear to lie on a straight line. 


What is the functional relationship between the size 
of the array and the number of draws of random numbers 
needed to fill that array? What values should be used 
to determine that function (e.g., are the standard 
deviations of any use)? How does one determine the 
type of function to be fitted? 


More data is needed, particularly for arrays larger 
than 20 x 20, 
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Standard 
Size Minimum Maximum Mean deviation 
4 28 9.27 4, 2897 
11 65 30.43 9.8916 
38 132 alee 18.4803 
69 2h9 132.73 32.7278 
120 397 231.41 54.3751 
230 622 355.39 67 4347 
334 874 517 42 105. 3863 
476 1201 718 .00 142.8735 
550 1405 906.27 159.409 


7T9T 1967 1192.07 208.446 

537 2230 1519.04 263.4189 
2885 1834.08 305.5707 
1512 3306 2287.07 340.8515 
1843 4221 2719.07 435.7215 
2172 4595 3276.58 528.3780 
2722 5206 3867.54 514.2681 
3229 6724 4498.46 §©6683..6533 
3727 7536 5255.13 797.8735 
4181 8697 6210.23 953.5309 
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For example, on a 10 x 10 array, all cells 

were hit in as few as 550 tries and at most 

1405 tries; the average number of tries was 7 
906 with a standard deviation of 159. 
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In our Contest No. 8 (PC39-4), each entry consists 
of nine 3-digit positive numbers in the range 100-999. 
To check the entries, a Fortran program is to be written 
to accept the nine numbers for one entry, punched on a 
card in columns 1-27. The program should reject, with 
an error message, any card for which the data numbers 
are not in the required range. 


The output from the program should show the input 
data and, on the same print line, the sum, the sum of 
squares, the mean (correct to 3 decimal places), and the 
sample standard deviation (correct to 3 decimal places). 
The End-of-file procedure is simply running out of input 
data cards. The following test data should produce the 
indicated results: 


Minimum FORTRAN 


Sum Sum of Mean Standard 
squares deviation 


100 101 102 103 104 105 106 107 108 936 97404 104.000 2.739 
100 200 300 400 500 600 700 800 900 4500 2850000 500.000 273.861 
500 501 502 099 504 505 506 507 508 UT MF RANGE 

500 501 502 503 504 505 506 507 508 4536 2286204 504.000 2.739 
123 234 345 456 567 678 789 890 901 4983 3404441 553.667 284.060 
999 999 999 999 999 100 100 100 100 5395 5030005 599.444 473.815 
123 123 123 123 123 123 123 123 999 1983 1119033 220.333 292.000@ 
123 456 789 987 654 321 234 345 456 4365 2732229 485.000 277.309 
991 992 993 994 995 996 997 998 999 8955 8910285 995.000 2.739 


All of the above is straightforward and is readily 
coded in Fortran. The Contest is this: using only ANSI 
Fortran, the $25 prize goes to the program using the 
least number of characters in the program, counting blanks. 
The number of characters per card will be counted up to 
the last card column used by the compiler. Thus, the 
Fortran statement: 


Contest If 
PROBLEM 138 


N=N+1 


would count as eleven characters, since the statement must 
begin in column 7 of the card. Similarly, the statement: 


45 IF(N.NE.NSAVE.@R.K.NE.KSAVE)GOTO10 


would count as 40 characters. 


The characters used on COMMENTS cards will not be 
counted. Each entry must include a listing of the progran, 
the test results, and the count of the characters. All 
entries must be received by POPULAR COMPUTING, Box 272, fe 
Calabasas, California 91302 by November 30, 1976. 
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In the May issue of 65 Notes, a program was given 
by Dr. Mordecai Schwartz to generate "Sequences." As 
implemented on an SR-52 calculator, the program will 
produce an endless string of sequences, of the sort 
found on intelligence test ("What are the next two terms 
of this sequence?"). Here are some examples of the 
output of Dr. Schwartz's ingenious program: 


Geeee. 57 155 “400. 1086 . 2801 7603 19608 
See e-o ehO ye —10 <-6 --20 -22  -46 <=64 

Pome oe cos e309 «32215 49193289 . 1159735 

Pee fot. 19-9 36° 53° 91. We. 235 
Temeetce 10G mea52 § 648" 1512 3888 9072 

2) 8 -8 -4 12 12 -28 -20 52 44 = -108 
Bpemeoe163 950° 5537 32272 188095 1096298 ... 


The program is refreshingly original and, even when 
you know exactly how it works, it still provides a 
fascinating game of solitaire. 


The flowchart on the following pages shows the logic 
of the program. 


Dr. Schwartz writes “The range of digits and 
their number is not accidental. The 12 integer- 
choices (-4 less than or equal to k less than or 
equal to 7) generated by RAND and utilized at 
various points in the program (see the flowchart) 
allow for the generation of interesting and varied 
sequences without a too-rapid increase in absolute 
values, for a bias in the direction of positive 
sequence terms, and for tight coding in the 9 
sequence selection process." 
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allow tight programming: 
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Contest 6 Results 


Contest No. 6 (PC37-4, April) called for the longest 
list of numbers which were multiples of consecutive integers 
and composed entirely of O's and l's. A sample list was 
given, with entries such as: 


18 times 61728395 = 1111111110 

19 5269 100111 

27 411522593 11111110011 

29 37969 1101101 

36 280558641975  10100111111100 


The contest was won hands down by Sam Wagner, Bluffton 
College, Bluffton, Ohio. Mr. Wagner chose to interpret the 
problem slightly differently, and produced a list of products 
of the first 1868 integers in which each product had an 
unbroken string of 1's followed by an unbroken string of O's. 


Thus, for the five numbers given in our example, Mr. Wagner 
shows: 


18 times 61728395 = 1111111110 
19 5847953216374269  1111213111111111111 
27 4115226337448559670781893 121211111111111 
1212111111111 
29 38314176245210727969348659  11111111111111 
112121111121111 
36 308641975 111211111100 


The listing of the 1868 results runs to 200 11 x 14 
pages, and includes items like: 


1864 times 59608965188 364 3299952312827849 3085 3600381497 
37720553171196948020982355746 304244158321411 
540295660467 3342870767 76 347 16261 325703385789 
22269909 394 372913686218407248450166905102527 
420123986647591797806 39008106819265617548879 
35145445875 

for a product of 232 ones followed by 3 zeros. 


In many cases, the second factor is 1820 digits long. 
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Comments on the article "ASSEMBLERS" in issue 41, page 8 
--by Irwin Greenwald 


At the risk of being offensive, let me begin by 
saying that the whole article strikes me as pomposity 
borne out of ignorance. Sorry, but you're simply out 
of date on the subject. Let me introduce you to the 
state of the art as of 1976. 


First, you are confused between the tasks of an 
assembler, the operating system it works under, the 
link editor, and the loader. Then, you do not seem to 
appreciate the trade-offs that the writer of an assembler 
must consider. For example, the choice of binding 
time is crucial to the design of an assembler. For each 
feature that you list as desirable, there is a price to 
be paid. Some desirable features will demand more 
storage during assembly, and storage may be in short 
supply; others would unduly lengthen the translation time. 


The writer of the assembler must seek some sort of 
optimum that will benefit most users, The point is that 
an assembler cannot be viewed in isolation; it interacts 
with many other programs, some of which are out of the 
control of the programmer of the assembler. The items 
that interact are: 


The source language. 

The (relocatable) object language. 
link editor. 

The loader. 

The operating system. 
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Each of these influences the others, and there are 
trade-offs among them. For example, decisions in the 
assembly process can be delayed if the object language 
(and the link editor) are powerful enough. Thus, one can 
write a one-pass assembler and relegate binding of forward 
references to the link editor provided that one insists 
that all external references are declared in the program 
prelude (at the front). But this is a restriction on the 
user. It buys him fast assembly, but a complicated (but 
not necessarily slow) link editor (it will still be I/0 
limited). The appropriate approach is obviously in the 
eye of the beholder. 
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A second aspect is the size (and number of passes) 
of the assembler, the speed of assembly, and the size of 
the source program it can handle. As the amount of main 
storage used for tables is increased, the size of source 
program that can be assembled is decreased (without 
resorting to table overlay). Hence, as the instruction 
dictionary is expanded, or the number of characters per 
symbol is increased, or the number of macro skeletons 
required per program, or the verbosity of error messages 
(assuming that they are kept in main storage), the size 
of the program that can be assembled is reduced. As we 
progress down the line from mini to micro, this issue 
becomes more important. 


The operating system is important in terms of the 
services it affords and restrictions it imposes on the 
assembler per se, as well as to the program being 
assembled. At a minimum, the assembler must provide 
linguistic elements (usually via a library of macros) 
that enable the programmer to call on the operating system 
services and specify exception handling (error and end 
conditions). Clearly, the syntax and semantics of these 
macros is dictated more by the operating system than by 
the writer of the assembler. 


I take special issue with your remark that no 
assembler has ever satisfied your list of desirable features. 
Most modern assemblers do even more (with the possible 
exception of verbose error messages) functionally than 
you have asked for. Most modern assemblers include the 
following features: 


1. The ability to specify your own mnemonic op-codes 
and synonyms for them, which is superior to any built-in 
synonym dictionary. 


2. Symbols with simple constraints: 


a. Less than N characters (N varies from 6 
on up). 


b. Restrictions on lead character (alpha or 
alphanumeric). 


ec. Restrictions on embedded characters, 
usually no blanks or arithmetic 
operators (+, -, *, /). 


3. Address arithmetic with, unfortunately, complex > 
rules for expressions involving relocatable 
symbols--a point which you have ignored. $ 


4, Detection of multiply defined, and undefined, 
@ symbols, with different error messages. 
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Plentiful pseudo-ops. 
Conversion of constants. 


Many list options. 


on nw 


Error messages, running the gamut from 
extremely concise to terribly verbose. 


9. Listings with line numbers keyed to the 
compressed deck order and update (ALTER) 
corrections based on the line number. 
(Continuity numbers are meaningless without 
SOS-style CHANGE corrections, which are almost 
impossible to design in a human-engineered way.) 


10. Commentary on instruction lines as well as 
independent commentary lines. 


11. Built-in (operating system) macros; macros (as 
well as code) callable from a library, using 
INCLUDE. 


@ 12. User-defined in-line, or library-called macros. 


13. Local symbols via various coding schemes (I 
would prefer different syntax and semantics 
in most cases). 


14. Zero level addressing, including specification 
of decimal, floating, hex or octal, and 
character constants. 
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4) 
is pointed out in Problem 7-11 of Litton's Problematical 
Recreations (edited by James Hurley, 1971). What other 
numbers A and B have the property that 
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Z-Sequence 
A-Sequence 
Me Sequence 


Dr. Mordecai Schwartz's work on Sequences (page 5 
of this issue) triggered research into other possible ways 
of generating interesting sequences. Dr. Schwartz's 
sequences are intended to be predictable; that is, given 
enough terms, one can readily deduce the pattern of 
generation. 


The twelve types of generator he used involve 
some mix of: 


...the two previous terms 
...the term number 
»..one or two constants 


but there is plenty of room within that framework to 
seek new sequences. Some of these turn out to be dull, 
such as: 


+a 


n n-1 n-2 =A 


However, the sequence: 
a, -|24,. +4). —A| 


turns out to be fascinating. On the facing page are the 
first 156 terms of the sequence (taking the first two 
starting terms each to be unity). The sequence (call it Z) 
has these properties: 
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28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
4 
42 
43 
44 
45 
46 
47 
48 
ug 
50 
51 
52 


DS 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
Te 
73 
74 
75 
76 
77 
78 


89 


105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 


131 
132 
133 
134 
135 
136 
137 
138 
139 
140 
141 
ye 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
155 
156 


PC42-13 


PC42-14 


1. Successive terms seem unpredictable; that is, 
the values Jump around, but not too wildly. 


2. In the long run, its terms tend toward n/3, 
where n is the term number, Thus, the sum of the first 
8002 terms is 11245744, whereas the sum of n/3 for 
n=l, 2, 3,...,8002 is 11213501, a difference of only 
32233. 


3. The sequence is readily generated in any 
language on any machine, using the simplest of operations, 
and dealing only with positive integers. The accompanying 
table gives restart values at various points in the sequence. 


4, No term after the first is ever as large as n. 
The circied numbers in the table show the appearances of 
new larger values. Such new larger values occur 68 times 
in the first 5081 terms. 


5. Zero appears at terms 3, 7, 11, 28, 31, 140, 
239, 600, and 6476, but does not seem to appear after that. 
Up to n = 5000, only 67 terms have a value less than 5. 


6. Successive equal terms appear at n = 90 and n = 98, 
but this phenomenon does not seem to occur again. 


7. If item 2 above holds indefinitely, then the sum 
of the series should be approximately the sum of n/3. 
The sum of each term divided by n should be approximately 
n/3 itself. The sum of each term divided by the square of 
the term number should be approximately 1/3 the sum of 
the !armonic series (see PC9-14). The sum of each term 
divided by the cube of the term_number, however, should 
approximate 1/3 the sum of (1/n@), which should converge. 
It appears to converge to 1.2346... 


8. Almost every integer appears as a term in the 
sequence sooner or later, except for a few numbers, such 
as 245, 449, 569, 575, and 903. 


Various statistics that have already 
built up about the Z-sequence are given 
on the following pages. Items 5, 6, 7, 
and 8 above indicate the need for further 
exploration of this new function. 


prosim 140 


Ld 


100 1530 2900 91450 5700 186112 
200 4976 3000 92000 5800 182604 
300 8286 3100 109134 5900 190098 
400 12220 3200 111310 6000 202632 
500 14702 3300 103512 6100 201228 
600 18510 3400 3 §=113304 6200 213454 
700 22532 3500 118204 6300 215678 
800 2hg72 3600 114396 6400 205668 
900 29744 3700 121572 6500 240600 
1000 33124 3800 122420 6600 209158 
1100 34840 3900 124572 6700 210156 
1200 37634 4o00)=3—-« 125240 6800 216194 
1300 43490 4100 =. 1.442410 6900 230826 
1400 47800 4200 »=—. 1142714 7000 230776 
1500 44gok 4300 140472 7100 208950 
1600 52124 4400 8=—. 149084 7200 86241294 
1700 54240 4500 =. 155144 7300 238912 
1800 60852 4600 8152844 7400 243070 
1900 59192 4700 §=©159518 7500 249480 
2000 66428 4800 156766 7600 256620 
2100 71988 4900 =. 162708 7700 270472 
2200 69482 5000 159268 7800 250886 
2300 77658 5100 177284 7900 8275186 
2400 82708 5200 174080 8000 267358 
2500 83586 5300 173936 8100 285022 
2600 89790 5400 184918 8200 262086 
2700 87008 5500 178608 8300 279974 
2800 88490 5600 177888 8400 268484 


Sums of successive groups of 100 terms of the 
Z-sequence. 


Each entry can be compared to the corresponding sum 
of (n/3) by adding the hundred fumes to the preceding 
hundred number and multiplying by (50/3). For example, 
the last entry compares to (akos sucom crane 275000. 
The discrepancy for the last entry is thus only 2.37%. 
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Table of restart values for the Z-sequence. The 
entry at 1500, for example, gives the values of the terms 
numbered 1498, 1499, and 1500. 


Thus, 1500 ~- 2(23) - 270 = 1184. 


The GILBREATH Conjecture 


The accompanying figure illustrates the Gilbreath 
conjecture.* The left-hand column contains the consecutive 
prime numbers. Subsequent columns contain the absolute 
values of the differences. The leading number in each 
column of differences is one, and Norman Gilbreath 
Conjectured that this phenomenon would hold good indefinitely. 
The article cited states that it has been verified for all 
primes less than 792 722 (that is, for the first 63 419 
primes). 


At first glance, it seems remarkable that the primes 
should behave so well, inasmuch as the primes distribute 
themselves irregularly among the natural numbers. What 
is really going on here begins with column 2, which has a 
leading 1 followed by small even numbers. For small even 
numbers, absolute differences will converge to 2 or 0 
rapidly as one goes across a difference table. The 
whole pattern hinges on having the entries in column 2 
relatively small in relation to what line they are on. 


Consider the last row of the figure. How large 
would the left-hand entry have to get to disturb the 1 in 
the far right position? Sticking to the primes, it could 
get to 83 and still satisfy the pattern. That is to say, 
an entry as large as 22 on the 16th line of the first 
differences would be tolerable, as long as the entries 
above it were smaller. The first occurrence of a 
difference of 22 in the primes is at 1129, 
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a difference of 6 23 first appears 
10 139 with the prime 
16 1831 shown here 
20 887 
30 4297 
40 19333 
50 31907 

100 396733 
150 13626257 
166 83751121 
200# 378043979 
220 47 326693 
3624 10726904659 


Source: "Statistics on the First Six Million Prime Numbers," 


F. Gruenberger and G. Armerding, October 1961, The RAND 


Corporation P-2460, for all entries except those marked (#), 


the latter due to T. R. Parkin and L. J. Lander, 1967. 


*In "On a Conjecture Concerning the Primes," by 
R. B. Killgrove and K. E. Ralston in Mathematical 
Tables and Other Aids to Computation, 1959, p. 121. 
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In other words, the Gilbreath conjecture simply 
underscores the fact that differences between successive 
primes grow, but very slowly. The table shows some of the 
first appearances of differences. Thus, a difference of 
100 will attenuate properly on the Gilbreath pattern, as 
long as there are about 100 columns of differences, which 
there would be at the line on which the prime 396733 
would appear. So the conjecture is not quite so 
remarkable, but still interesting. 


Which leads to a problem. Suppose we fabricate a 
new pattern as follows. Generate for the second column 
the sequence shown at the right: 


This is exceptional:— 
and we extend this sequence to 1000 lines 


(the last entry will be 18, in the set 
beginning 2, 4, 6, 8, at line 992). 


ae eee eee 


This pattern will satisfy the conjecture 
also; that is, if extended for absolute differences, 
the leading number for every column of differences 
will be one. 


The problem is: how large can that last 
entry be before the pattern is spoiled? 
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Log 42 1.62324929039790046 322098 305657 2244529451891141976770 
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An 42 3.737669618283368305917830101823882002 36007 5421764866 
V2 6.4807406984078602 309659674 36087996657 705204 307058347 
2 3.4760266448864497867 39865219004 5374 340048385 38786967 
VE2 1.453198460282267816555796587579 32369821041 37217 30762 
'VE2 1.03808398947 16 3267547 312988 37 207913941284 25076802692 


e 173927494.1520501047 . 394681303611235226147984057725008 
40103706101184 3950097817 387725921 


7 759092417 205229 39087 3 .876361772680695986 3134896756244 
89165791115324 37010794474495291 


tan-1 42 1.54699130060982667577745 36 36281960677419 304556603822 
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| Sf anything ran go wrong, it will 


...and it will happen to you 


-.. and af the worst possible time 
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